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Abstract 

We present kinematical distributions for top quark pairs produced at the 
Tevatron pp Collider, as predicted within Next-to-Leading-Order QCD. We 
consider single and double-inclusive distributions, and compare our results 
to those obtained with the shower Monte Carlo HERWIG. We discuss the 
implications of our findings for experimental issues such as the measurement 
of the top quark mass. 
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Now that the existence of the top quark has been firmly established via its 
detection in hadronic collisions |], @, |3|, experimental studies will focus on the 
determination of its properties. In particular, the measurement of its mass and of the 
production cross section and distributions will certainly be among the first studies 
of interest. The top mass represents one of the crucial parameters for testing the 
Standard Model ||, while the production properties, should they display anomalies, 
could point to the existence of exotic phenomena |J. 

Early estimates of the total production cross section for top quarks || 0] have 
later been updated ||. In ref. || a partial resummation of leading logarithms from 
soft gluon emission has also been included. 

In this Letter we present a set of kinematical distributions of relevance for the 
study of production dynamics, as obtained from full Next-to- Leading-Order (NLO) 
QCD perturbative calculations. We will consider inclusive p T and pseudorapidity 
(if) distributions, as well as invariant mass (M(tt)), transverse momentum (Pf ) 
and azimuthal correlations (A0) of the top quark pair. 

The inclusive p T distribution is sensitive to channels such as Wg-^tb [[II]], which 
are found to contribute with a small cross section, predominantly at low p T . The 
invariant mass of the pair is an obvious probe of the possible existence of strongly 
coupled exotic resonances, such as technimesons 0. The transverse momentum 
of the pair is an indication of the emission of hard hadronic jets in addition to 
those coming from the top decays. The presence of these additional jets generates 
potentially large combinatorial backgrounds to the reconstruction of the top mass 
peak from the decay products . An accurate understanding of these backgrounds 
is in principle very important for a precise measurement of the top mass. 

The NLO total, single and double inclusive cross sections for the production 
of heavy quarks in hadronic collisions have been calculated in references ||, JTT[ 



and fll2"l , respectively. Since most experimental studies are performed using shower 



Monte Carlo event generators |13| to simulate the behaviour of the events in the 
detectors, we will compare our NLO results with those we obtained with HERWIG 
|14|| . This is important in order to assess the reliability of the theoretical inputs used 



by the experiments. It has been recently pointed out |15| that there are possibly 



significant differences in the jet activity predicted by fixed-order parton level calcu- 
lations and the results of HERWIG, with possible implications for the experimental 
mass measurement. We will discuss here a possible origin of these discrepancies. 

For our calculations we use m t =176 GeV (CDF, ref. [1]]). Differences of ± 
15 GeV, approximately the current experimental error, do not affect the shapes of 
the distributions we present nor any of our conclusions. Our default set of parton 



densities is MRS A [16|, and we explore a range of renormalization and factorization 
scales /j /2 < /i R = /j f < 2/i , with /i 2 = m t 2 + p T 2 (/V = m t 2 + {{p T ) 2 t + (p T ) 2 )/2) 
for single (double) differential distributions. 

Throughout our plots, we rescale the HERWIG calculations by the perturbative 
K factor given by the ratio of the 0(a^) and the 0(a 2 ) results. The K factor is of 
the order of 1.3 for all choices of parameters. The distributions of p T , r\ and M(tt) 
are shown in figures 1-3. The solid lines correspond to the NLO results obtained 
using /i R =/i F =/i - The curves obtained using a different choice of scales, within the 



1 



10" 



10 



ft 



ft 

13 



b 



10 



10 



10" 



-2 



-3 _ 



-4 _ 




Solid: NLO QCD 
o : K * HERWIG 

,2 i „2 



/i=/x = V(mj+pT), MRSA 
m t =176 GeV 



100 



200 



300 



Pt (GeV) 



Figure 1: Inclusive p T distribution of the top quark. 



range quoted earlier, result in an overall normalization change, by approximately 
±10%. No change in the shapes is observed, and we therefore decided not to include 
these curves in the plots for the sake of clarity. The square points correspond to the 
HERWIG prediction for /i R =/i , rescaled by a K factor equal to 1.34. 

The figures indicate clearly that distributions which are not trivial at leading 
order, such as p T , r\ and M(ti), agree perfectly between the NLO and HERWIG 
calculations, at least for values of p T and M(tt) not too large. This result follows 
naturally both from the observed stability of those distributions under radiative 



corrections regardless of the heavy quark mass [12| , and from the large inertia of the 
top quark against the combination of higher order and non-perturbative corrections 
introduced by the Monte Carlo. At large p T and M(tt), however, one starts observing 
a small deviation, due to the multiple gluon emission from the final state top quarks. 
This effect is not described by the NLO QCD calculation, which only includes one 
gluon emission. Multiple gluon emission becomes relevant for p T ^> m t , where it 
leads to a top quark fragmentation function softer than predicted by the fixed order 
NLO calculation. Due to the smallness of the cross section in these regions, large 
statistics will be needed in order to test this behaviour. 

Those distributions which are trivial at leading order, A0 and P\ , are viceversa 
most sensitive to multiple gluon emission from the initial state. This is because even 
small perturbations can smear a distribution which at leading order is represented 

by a delta function, which is the case of and A0. The most visible effect is 

observed in P~ , Figure 4, where we include the NLO curves relative to the three 
choices of scales, /i H =yU F =/io (solid), /xo/2 (dots) and 2// (dashes). Contrary to the 
previous cases, significant differences in shape arise here among the three choices in 
the small p T region. The HERWIG result (normalized to the area of the solid curve) 
is also shown. The NLO and the HERWIG distributions assume the same shape 
only for P,jr larger than approximately 20 GeV. 
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Figure 2: Inclusive pseudorapidity distribution of the top quark. 





(P**) (GeV) 


F(P** > 10 GeV) (%) 


F(P** > 20 GeV) (%) 


NLO QCD, /i R =/i F =/i 


11.9 


30 


15 


NLO QCD, ji R =f/ P =2^ 


9.1 


23 


11 


NLO QCD, n K =fji F =fjLo/2 


16.6 


44 


22 


HERWIG, /i R =/i F =/i 


17.5 


51 


28 



Table 1: Average value of P^ and fraction of the cross section with P^ 1 > 10 and 
20 GeV for different calculations. 



To better quantify the differences, we quote in Table 1 the value of the average 
and of the fraction of the total cross section with P^ > 10, 20 GeV for the 
different curves appearing on the plot. First of all we notice the strong scale 
dependence of the results in Table 1. This is a first indication of the lack of reliability 
of the fixed order calculation for these particular quantities. Observe that the NLO 
result obtained with the lowest scale is closer to the HERWIG result. In fact, a lower 
scale is presumably more appropriate in this problem. This is because there are two 
relevant scales in this case: the mass of the top, and the transverse momentum 
cut. Since the production of top quarks takes place mainly through the exchange 
of a highly virtual gluon (with virtuality of the order of the top mass), while the 
emission of the jet takes place at the smaller scale of the jet transverse momentum, 
it is sensible to take in the 0{a\) cross section two powers of a s at the large scale /io 
and one at the transverse momentum of the jet. We have performed this exercise, 
and obtained a(P^ > 10 GeV) = 1.987 pb, which is 43% of the total cross section 
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Figure 3: Invariant mass distribution of the tt pair. 



computed at the scale fi . 

The high value of (P? ) for the NLO calculation indicates that the probability 
of emitting a gluon with transverse energies k T of the order of 10 GeV is of or- 
der 1. Technically, this is related to the appearance of a large double logarithm 
0(a s log 2 [m 2 //c^] ) . It is well known |I7J that in these circumstances multiple gluon 
emission, corresponding to multiple powers of these double logs, needs to be re- 
summed. A typical example can be found in W and Z production at small k 7 



18 



The Monte Carlo HERWIG performs the resummation of leading and some sub- 
leading double logarithms, faithfully reproducing the behaviour of the analytically 
resummed cross sections in known cases |14[] . The large top mass extends the region 
where Sudakov effects are important up to k T values which are often considered 
within the domain of applicability of standard fixed-order perturbation theory. 

We have implemented a leading logarithm version of the resummation procedure 
described by Collins Soper and Sterman in |18|, and qualitatively reproduced the 



behaviour of the Monte Carlcf]. We therefore conclude that the parton level result 
(with the factorization and renormalization scale taken of the order of the mass of 

the top quark) is less reliable than the Monte Carlo in the small P^ region, and 

gives a cross section (for P^> 10 GeV) which is roughly 60% of what one would 
get with a more judicious choice of the scale, and with the inclusion of Sudakov 
effects. We also point out that no choice of scales in the NLO calculation would 
describe the detailed shape in the region P| < 20 GeV. On the other hand, we 

expect the parton calculation to be more accurate in the large region, where the 
soft approximation used in the Monte Carlo breaks down. 

A similar problem should be expected in the description of jet radiation during 
the production of any massive system, in particular W or Z boson pairs. This might 
affect parton level calculations of backgrounds to tt production in the dilepton plus 



3 A similar calculation in the context of b pair production has been given in rcf. |TJ 
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Figure 4: Transverse momentum distribution of the tt pair. 



two jet channel. A simple estimate along the lines of what performed above for 
the top, sets the threshold below which fixed order perturbation theory becomes 
unreliable at about 10-15 GeV, which still corresponds to jets sufficiently energetic 
to be reconstructed as such by the Collider detectors. A confirmation of this estimate 
can be found in the study of resummation effects on the p T distribution of Z pairs 
presented in reference |20| . 

The observations made above partially explain the disagreement found in ref. 
|T5|] between the relative fraction of jets emitted during the production process and 



jets emitted from the b quark resulting from the top decay. These authors found 
that their parton level calculation gave rise to a ratio smaller than HERWIG's. They 
interpreted it as a deficit in jet radiation off the final state b quark within HERWIG. 



In order to make a more definite comparison with the results of ref. [|T^], we have 
computed the ti jet cross section with the same value of the parameters and the 
same set of cuts that they applied. We used therefore m t = 174 GeV, /io = rn t , 
\r)j e t\ < 2.5, > 10 GeV. In order to emulate the cuts they perform on the b jets 
coming from top decays, we let the top quark decay isotropically, and imposed the 
additional cuts [ 77^ | < 2.5 and > 10 GeV. We verified that these additional cuts 
affect the cross section by less than 3%. We obtained a cross section of 1.30 pb, to 
be compared with the value of 1.0 pb of ref. [T5|]. We have no explanation for this 
difference. From the considerations made previously, we expect that the inclusion 
of Sudakov effects should increase the cross section by roughly 66%, leading to a 



value of 2.16 pb, more than a factor of two larger than the result of ref. |T5 . 

Figure 5, finally, shows the tt azimuthal correlations. Notice that most of the 
shape discrepancies observed at P^<20 GeV are forced here into the bin around 
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Figure 5: Azimuthal distance distribution of the tt pair. 



the A0=7r point. This is because even a 20 GeV kick will only deflect the trajectory 
of the t or t quark (each traveling with momentum of the order of m t /2), by about 
1/10 of a radian. 

While completing this paper, a study of single inclusive distributions of top 
quarks at the Tevatron appeared |HJ] . These authors improved the NLO calculation 
with the inclusion of a partial resummation of leading soft logarithms ||. Their 
results indicate that single inclusive distributions are only affected in overall nor- 
malization, but not in shape (at least for p T not much larger than m t ). This is 
consistent with the results presented here, and provides yet more confidence in the 
validity of the perturbative predictions applied to single inclusive quantities. 

In concusion, we presented here a set of kinematic quantities which should be 
used to test the dynamics of top quark production in hadronic collisions. Most 
of these quantities were shown to be subject to very small theoretical uncertainty 
within NLO perturbative QCD, and to be rather insensitive to higher order or 
hadronization effects. One noticeable exception is the transverse momentum of 
the pair, which is in turn related to the spectrum of jets recoiling against the top 
pair. We found in our study that an accurate description of the region P**< 20 
GeV requires the resummation of leading soft and collinear double logarithms, as 
implemented in appropriate shower Monte Carlo programs. 
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